Spin current through an ESR quantum dot: A real time study 
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The spin transport in a strongly interacting spin-pump nano-device is studied using the time- 
dependent variational-matrix-product-state (VMPS) approach. The precession magnetic field gen- 
erates a dissipationless spin current through the quantum dot. We compute the real time spin 
current away from the equilibrium condition. Both transient and stationary states are reached in 
the simulation. The essentially exact results are compared with those from the Hartree-Fock ap- 
proximation (HFA). It is found that correlation effect on the physical quantities at quasi-steady 
state are captured well by the HFA for small interaction strength. However the HFA misses many 
features in the real time dynamics. Results reported here may shed light on the understanding of 
the ultra-fast processes as well as the interplay of the non-equilibrium and strongly correlated effect 
in the transport properties. 
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I. INTRODUCTION 

The rapid progress of nano-electronics and information 
technologies has prompted intense interest in exploiting 
the spin properties of the electrons, which results in the 
emergence of spintronics.[lJ One of the most important 
spin-based electronic devices is a mesoscopic quantum 
dot (QD) system. Spin-polarized transport through a 
QD has been extensively investigated recently. It has 
been shown theoretically [2] and demonstrated experi- 
mentally [3] that a QD system will function as a phase- 
coherent spin pump in the presence of sizable Zeeman 
splitting. Very recently, spin-polarized current has been 
detected from a quantum point contact (QPC) [4] and 
from Coulomb-blocked QDs |5J. The electron spin res- 
onance QD (ESR-QD) could also serve as an element 
device for quantum computing |6J. 

Due to its small size, Coulomb correlation could play 
important role in the transport experiments involving a 
QD. At low temperature, Kondo effect creates new states 
with many-body character at the Fermi level. Although 
the effect of Kondo resonance on the charge current has 
been studied in the QD for different situations, its in- 
fluence on the spin current is less known. The physi- 
cal processes involved in the transport experiments are 
out of equilibrium in many cases. Moreover, the state of 
the system may also be time dependent. These features 
make the investigations unaccessible from the conven- 
tional many-body tools. So an essentially exact numeri- 
cal method for the non-equilibrium and time dependent 
phenomena in the interacting nano-devices is highly de- 
sirable, which can also verify the approximations used in 
various analytical approaches. 

There exist several powerful methods to deal with the 
low dimensional correlated systems, such as the numeri- 
cal renormalization group (NRG) and the density-matrix 
renormalization group (DMRG). With the input from 
the quantum information science, time-evolving block 
decimation (TEBD) [8] and adaptive time dependent 



DMRG (t-DMRG) have received much attention. \W\ 
To our knowledge, there are some previous studies of the 
non-equilibrium transport of the nano-devices using the 
adaptive t-DMRG technique. [TTHT3] By adopting the log- 
arithmic discretization, da Silva et al. [14J have studied 
the Kondo correlations. Guo et al. [15] have examined the 
noninteracting resonant level model in the Landau- Zener 
potentials. Besides the adaptive tDMRG approach, there 
are also attempts based on the time-dependent NRG 
|16| , functional RG [17], Dyson equation embedding |18| , 
flow-equations [19, 20J, and quantum Monte Carlo [2TH23] 
methods. But the adaptive t-DMRG approach gives di- 
rect access to the transient regime and could handle the 
time dependent Hamiltonian directly. The tVMPS ap- 
proach adopted in this paper is directly connected with 
the adaptive t-DMRG [TO] and the TEBD [3 |8] ap- 
proaches. The computational cost of these two methods 
is very similar, and in practice both methods achieve a 
similar accuracy. [24J. 

The merit of the VMPS approach lies in two aspects. 
First, it represents a large class of states, which could be 
seen from the success of the NRG and DMRG approaches 
(which generate MPS in their processes) for the zero and 
one-dimensional quantum models. And the fact that the 
entanglement entropy increases slowly in one dimension 
also permits one to simulate the states classically using 
the VMPS method. [24J. Second, it is easier to handle 
the MPS i.e., the overlap of two MPS, the expected value 
of an operator in a given MPS etc. can be calculated with 
polynomial complexity. 

In this paper we study the spin current through an 
ESR quantum dot with Coulomb interaction [25H27] . 
We obtain the transient as well as the quasi-steady 
state spin current using the time-dependent VMPS 
method. We also studied the effect of the interact- 
ing on the spin current. The results are compared 
with time-dependent Hatree-Fock approximation (TD- 
HF) and non-equilibrium Green's function (NEGF) ap- 
proach for the quasi- steady state. 
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Figure 1: (a). The experimental set up of an ESR quantum 
dot. Besides the Zeeman field Bo along the z-axis, there is a 
rotational magnetic field B\ acting in the x?/-plane. The dot is 
coupled to a noninteracting lead and forms a single impurity 
Anderson model (SIAM) chain. The SIAM chain is unfolded 
into two spinless chains, and they are coupled at the leftmost 
end by Coulomb interaction U and the rotational magnetic 
field B\. (b). Schematic view of time independent energy 
levels in the rotational reference frame. The Fermi level of 
the spin up/down electrons in the lead is shifted to zbcj/2, 
while the | «— ) and | —>) levels in the dot split by 2g\iBB\. 
The dot levels could be tuned by the gate voltage V g . 



o 8site, up 
o ibsite, up 




2 4 6 8 10 
t [Nt] 

Figure 2: (Color online) The spin current for different chain 
length L = 8, 16,32. The red solid line is the tVMPS result 
and the dots are from direct integration of the von-Neumann 
equation, see text. Quasi-steady state could be recognized 
for larger chain length. Where (hollow dots) and J± (filled 
dots) are of opposite sign. 



II. MODEL 

We consider an ESR setup of a quantum dot, where 
single-electron level of the dot is split by the Zeeman 
field Bq and the two spin levels are coupled by a rotating 
magnetic field £?i(cos(u;t), sm(ur)). The Hamiltonian 

reads H = H^ot + ^rotate + #lead 



Hdot = V g ^2 n ° Y~ ~ (n^-n^ + U n t n ; 

a 

^rotate = - 9 J^(d\d^ ^ d\d^~^) 
Hlead = -t'y^(<4- C l,cr + h.C.) 

cr 

N lead 

~ l (4,c^+l,a + h.C.) . (1) 

i=l;a 

Here d and c denote the annihilation operator of elec- 
tron in the dot and the lead. i^i ea d contains the terms 
describing coupling of the dot with the lead and the hop- 
ping in the lead. H rotate contains the rotating magnetic 
field, i^dot contains the Zeeman splitting, the gate volt- 
age terms and the on site Coulomb repulsion between the 
spin up and down elections. In the present study we fix 
t = 1 and gfiB^o = u = 1 where the ESR resonance con- 
dition is satisfied. We set gpB^i = 2 and t f =0A unless 
mentioned. The single lead is a noninteracting chain with 
Niead sites. The coupling of it with the QD is described 
by the hybridization T = nt' 2 p. 



III. METHODS 

The time dependence of the Hamilton could be 
eliminated by the unitary transformation U = 

e -iW2)E i « i c 4li +ct t c 4 , t )+(«|t di _4 dt )]_ ^ transforms 

the Hamiltonian into rotating reference frame (RF), see 
Fig[TJb). One can see that the rotating magnetic field, in 
effect, shifts lead electron energy to the opposite direc- 
tions for up and down spins. Quasi- steady state spin cur- 
rent has been studied in the RF using the non-equilibrium 
Green's function (NEGF) method in the noninteracting 
and infinite U case [25] and by NRG approach [27j in adi- 
abatic limit. In the following we will revisit the problem 
in transient as well as quasi-steady regime, and treat the 
interacting non-perturbatively with the exact numerical 
methods. 

For the time dependent Hamiltonian H(r), the evo- 
lution of the density matrix p follows the von-Neumann 
equation = [H,p]. For the noninteracting case, one 
could write H(r) with the single particle basis. Then 
direct integration of the von-Neumann equation with a 
given initial condition p(r)| r= o = po could be easily done. 
The initial density matrix po is calculated form ground 
state of the SIAM chain without the rotating magnetic 
field Si. 

From the density operator p{r) the occupation on 
the dot is calculated by n a (r) = (^(r^d^d^ir)) , 
and the spin current through the dot is J a = 
^Im[t / ('0(r)|(iJ.ci 5Cr |'0(r))], i.e. they are evaluated on the 
bond connect the dot and the lead. We set gpB = e = 
h = 1 in this paper. 

For the interacting case we adopt the VMPS approach 
to calculate the ground state of H(r = 0). Then we apply 
the rotational magnetic filed B\ on the dot at r = 0, and 



3 




O 5 10 15 20 25 30 35 40 

T [fl/t] 



Figure 3: (color online) Dot occupancy for each spin (upper 
panel) and spin current (lower panel) for a 64 sites chain with 
different gate voltage. The shaded region is averaged to cal- 
culate the quasi-steady state spin current. 



perform the time evolution. [28J. To reduce the dimension 
of the local Hilbert space, an unfolded technique [29] is 
used, the original SIAM chain is unfolded into two chains 
with different spins. The total length of the unfolded 
SIAM chain is L. They are connected at the end point 
by the Coulomb repulsion and the rotating magnetic field 
B\ . The errors of the computation mainly come from the 
following sources. First the Trotter decomposition error. 
Second the truncation errors accumulated in the course 
of the time evolution. For short time scale the Trotter 
error dominates while for the long time the truncation 
error dominates. Caution must also be taken because of 
the finite size of the leads. The electrons may bounce at 
the end point and the spin current flows along the reverse 
direction. This is an artifact of the present method and 
could be eliminated by careful finite size scaling analysis. 

We also use time-dependent Hartree-Fock (TD-HF) 
method to investigate the interacting case approxi- 
mately, in which the interaction term is factorized into 
U((n^)n^ -\-n^(n_i)). At each time step of the integration 
of the von-Neumann equation, the dot occupation which 
is used to update the Hartree-Fock Hamiltonian. This 
approach is as efficient as the non-interacting cases. 



IV. RESULTS 

In Fig. [2] we show the development of the spin current 
after the rotational magnetic field Bi is applied to the 
noninteracting dot. To reduce the transient region, we 
use a larger coupling constant t f = 1 here. For small chain 
length (L = 8, 16) the finite size effect shows up within 
the maximum time of our simulation. The charge pulse 
reaches the end of the chain and bounces back thus causes 



the reversing of the sign of the spin current. However, the 
spin current reaches a saturation value for larger chain 
length (L = 8, 16) which indicates the quasi- steady state 
spin transport is achieved. In the quasi-steady state, the 
spin up electrons flow in the dot, flip their spins and 
then flow away. There is no net charge transport since 
J-f = —J±. Overshooting behavior at short time scales 
(r ~ 1.8) is observed. It is due to abrupt applying of the 
rotational magnetic field and could be suppressed if the 
rotational magnetic field is turned on adiabatically. The 
tVMPS result for L — 32 nicely follows the direct inte- 
gration of the von-Neumann equation. The oscillation at 
long times in the tVMPS approach was also noticed in 
the previous study with the adaptive t-DMRG approach 
|13| . With the increasing of the bond dimension, the 
oscillation tends to disappear. The coincidence validates 
the VMPS method for the noninteracting case. However, 
its main power lies in the interacting cases and we expect 
similar precision could be achieved in that case. 

Electron occupation and spin current through the 
quantum dot are calculated as functions of time for dif- 
ferent gate voltages V g , shown in Fig 131 It is seen that 
the overall occupation on the dot oscillates as it reaches 
its quasi- steady- state value, the spin current develops in 
the mean time. Note that the transient state current 
can even be of the opposite sign with its value in the 
quasi-steady state, Figj3^b). By taking average of physi- 
cal quantities in the quasi-steady state, we extract vari- 
ation of the quasi-steady-state dot occupation and spin 
current with the gate voltage V g , see FigjEJ The spin po- 
larization in the steady state is less pronounced than the 
initial state. And it is even inverted for — 1 < V g < 1, 
Figj5|a). The steady state current attains its maximum 
value at V g = ±1. This is the case where = n\i.e. 
the dot is spin unpolarized. It can be seen from FigJTJb), 
in this case V g d= B\ = and the | <—) (| — »•)) is tuned to 
zero energy point. Spin flip process on the dot is most 
efficient and gives a maximum spin current. We also use 
NEGF to study the current and occupation in the RF, 
see inset of Figj5|a). It gives qualitatively similar re- 
sults, the small discrepancy is due to the finite size effect 
of the tVMPS approach. It should be noted that NEGF 
approach adopted here does not yield results of transient 
state since the unitary transformation eliminates the time 
dependence. To take into account the initial condition 
properly, the two-time Green's function should be used 
[30H32]. 

For interacting case, we calculate the real time spin 
current through the ESR-QD by tVMPS method. The 
results were compared with TD-HF results, see Fig(4]for 
U = 1. Although there are discrepancies in the real time 
data, the TD-HF approximation captures the overall be- 
havior of the spin current. This is due to the Kondo 
physics-which is missing from the TD-HF approach-does 
not manifests itself here since the presence of large mag- 
netic field B\ and Fermi surface splitting in the lead. 
Kondo effect may restored when B\ is reduced below 
Kondo temperature [27] . But small B\ also reduce the 
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Figure 4: (color online) The comparison of the tVMPS (full 
line) and TD-HF (dot line) results of spin current for different 
gate voltage, where the interacting strength U = 1. The 
shaded area is averaged to give the quasi-steady state current. 



chance of spin flip process on the QD, thus reduce the 
spin current through it. Surprisingly, the average of phys- 
ical quantities in the quasi-steady state of these two ap- 
proaches are in good agreement, see Fig(5j Based on 
this investigation, we validate that the HFA is a good 
approximation for the qualitative investigation of spin- 
pump devises away from the Kondo regime. However, 
caution must be taken when it is used to make predic- 
tion on the real time dynamics. For example for the 
V g = —1.0 case, the transient state spin current of TD- 
HF is of the opposite sign to the tVMPS prediction, thus 
it is quantitative wrong. 

The magnitude of spin current with the gate voltage 
is still a two-peaked curve similar to the noninteracting 
case. Interaction has different effects on the high and 
low filling regimes of the dot, Fi^5](b). It does not mod- 
ify the dot occupation and the spin current dramatically 
for V g > 0. The peak of the spin current remains at 
V g = 1 and the maximum value is only suppressed by 
one percent up to U = 3. However, since the average dot 
occupation is larger than 1 for V g < 0, the correlation 
effect shifts the spin current peak downwards by U. This 
is a manifestation of the Coulomb blockade effect. 

The fluctuation of the density A = (Sn^-Sni) = 
(n^n^) — (n^)(n±) is also calculated. It is a measure of the 
accuracy of the HF approach. This quantity is conserved 
under the particle- hole transformation. Thus two differ- 
ent gate voltages satisfying V g + Vg = — U give the same 
value of A. This is respected in the tVMPS calculation. 
But the discrepancy from the tVMPS is unrelated to it, 
i.e., although V g = 1 and V g = has the same A, the 
discrepancy of the spin current from the tVMPS result is 
not the same. We find the HF approximation is surpris- 
ingly good for the quasi-steady state physical quantities 



Figure 5: (Color online) Occupation of electrons on the dot 
(upper panel) and spin current through the dot (lower panel) 
for different gate voltages. The spin current of high filling 
region is not affected by the interacting, while the spin cur- 
rent peak at low filling is shifted by U. Inset are the U — 
quantities calculated in rotational reference frame by NEGF 
approach. 



even if the fluctuation is relatively large. The fluctuation 
is always negative due to the Cauchy-Schwarz inequal- 
ity, indicating that the HFA overestimates the potential 
energy. The absolute value of the fluctuation reaches its 
maximum 0.2 for — 1 < V g < 0. However, even for these 
cases, the HF data still shows good agreement with the 
tVMPS data. 



V. CONCLUSION 

To conclude, we perform essentially exact real time 
calculation of the spin current through an interacting 
ESR quantum dot. We benchmark the essentially ex- 
act tDMRG result against those from various analyti- 
cal and approximate methods. From the extracted aver- 
age spin current, we obtain the Coulomb block shift of 
the spin current peak, confirmed by the time-dependent 
Hartree-Fock calculations. We find that the spin current 
attains its maximum for a spin neutral quantum dot. The 
spin neutral condition is fulfilled for two gate voltages 
where V g matches the magnitude of the rotational mag- 
netic field. The two spin current peaks respond differ- 
ently to the electronic correlation, the lower filling peak 
shifts downwards by U while the higher filling peak is 
nearly unaffected. Comparison to the NEGF approach 
with an infinite lead shows that the finite size effect does 
not affect the qualitative behavior of the quasi-steady- 
state quantities. These results are also compared with 
those of the TD-HF approach. It is shown that the TD- 
HF gives accurate quasi-steady-state dot occupation and 
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spin current. However its prediction on the real time 
dynamics is problematic. 
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